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We study the influence of long-range interatomic interactions on the properties of supersonic 
pulse solitons in anharmonic chains. We show that in the case of ultra-long-range (e.g., screened 
Coulomb) interactions three different types of pulse solitons coexist in a certain velocity interval: 
one type is unstable but the two others are stable. The high-energy stable soliton is broad and can 
be described in the quasicontinuum approximation. But the low-energy stable soliton consists of 
two components, short-range and long-range ones, and can be considered as a bound state of these 
components. 
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As is well known W,0| , anharmonic chains with interac- 
tions between nearest neighbors can bear pulse solitons, 
compressive localized excitations which are very robust 
and propagate with supersonic velocities without energy 
loss. Because of their coherence, the solitons play an im- 
portant role in determination of dynamical, thermody- 
namic, and transport properties of one-dimensional an- 
harmonic systems ||l[. Among other things, they have 
been invoked in order to explain energy transport in DNA 

i- 

However, the interatomic interactions in real systems 
are strictly speaking long-ranged. In particular, the DNA 
molecule contains charged groups with Coulomb inter- 
actions between them |j]. Therefore, it is essential to 
clarify how the long-range interactions (LRFs) can affect 
the soliton features. It is generally believed that such in- 
teractions are very small (in comparison with the anhar- 
monic interactions between nearest neighbors) and can 
be safely neglected. However, as we show in the present 
paper, even very weak LRI's cause new qualitative ef- 
fects if the interactions are ultra-long-ranged. A striking 
illustration is a chain with pure (not screened) Coulomb 
interactions between charged particles where the sound 
velocity is infinite regardless of the intensity of these in- 
teractions. As a consequence the pulse solitons merely 
do not exist in such a model (whereas the pure Coulomb 
interactions do not prevent M the existence of immobile 
intrinsic localized modes therein). Generally, arbitrary 
LRI's introduce into the system a new length scale, the 
so-called radius of the LRI's. If the radius of the LRI's 
far exceeds the interatomic distance, the competition be- 
tween the length scales manifests itself in a number of 
qualitative effects (see Refs. |6|j^ for the exponential-law 
LRI's and Refs. ||-|l^ for the power-law LRI's). The 
greater is the radius of the LRI's, the more pronounced 
are these effects. 

In this paper we show that two types of stable pulse 



solitons can coexist in a certain interval of velocities in 
anharmonic chains with ultra-long-range interatomic in- 
teractions even if they are very weak. 

Let us consider a chain of equally spaced particles 
of unit mass whose displacements from equilibrium are 
Unit) and the equilibrium spacings are unity. The Hamil- 
tonian of the system is given by 

+ 2 X! "^ryi.riiUm " Unf > , (1) 

m'>n ^ 

with the anharmonic interactions V(w) = vp' j^ ~ 
w^ /?i between nearest neighbors and the harmonic LRI's 
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|to — n|^ between all particles 



of the chain. Here J characterizes the intensity of the 
LRI's whereas a and s determine their inverse radius. 
The parameters a and s are introduced to cover differ- 
ent physical situations from the limit of nearest-neighbor 
interactions (a 2> 1 or s ^ 1) to the limit of ultra- long- 
range interactions (a <C 1 and s < 3). The Hamiltonian 
(0) generates equations of motion of the form 
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w„ are relative displacements and 
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where w„ = Un+ 

F{w) = dV{w)/dw — w — w 

We assume in what follows that a 7^ (the case a — 
has already been considered in Ref. H); in doing so we 
studied most extensively two cases: the physically im- 
portant screened Coulomb interactions (s = 3) and the 
Kac-Baker LRI's (s = 0). However, in view of the fact 
(tested numerically) that all cases with < s < 3 lead to 
qualitatively the same results but the case s = allows 



also analytical consideration, we discuss only the case 
s = from this point on. 

In the quasicontinuum limit, treating n as a contin- 
uous variable \n ^ x , Wn(t) — > w(x,t), Wm(t) -^ 



;(x,t)] and keeping formally all terms in the 



g(m-n)9^ 

Taylor expansion of e^™^")^^ , the equation of motion 
for s = can be cast in the operator form 

[dl - JQ{a,d^)]w{x,t) - Asinh" (^]f{w) = , (3) 



where 



Q(a,9,) = (e" + 1) 



4sinh^(9:r/2) 



(4) 



with K = 2 sinh(a/2) is a linear pseudo-differential op- 
erator. The speed of sound c (which is an upper limit 
of the group velocity of linear waves), determined by the 
expression c^ = 1 + J(l -I- e^")/(l — e~")^ , grows indef- 
initely as a decreases. 

We are interested in the stationary soliton solutions 
w(x,t) = mix — vt) propagating with velocity v. In 
this way we reduce our problem to a nonlinear eigen- 
value problem with v being a spectral parameter. In- 
deed, substituting z = x — vt and using the continuum 
approximation 4sinh^(9a;/2) w 9^, we can write Eq. (||) 
in the form 



{dl - sl){dl sl)w{z) = ^^{dl - n^)w\z) 
where the parameters s± are given by 



(5) 




(6) 



The parameter s^ is finite at all velocities v > c and 
tends to \/l2 for v — > oo. The parameter s_ vanishes 
at V = c and tends to k for w — > oo. Using the Green's 
function method ||9|] one can show that stationary soliton 
solutions exist only for supersonic velocities v > c. The 
properties of these solitons are determined by the ratio 
of s+ and s_. In Fig. H we plot the energy of the soliton 
solutions of Eq. (|^) which were found numerically using 
the method developed in Ref. |^. 

The soliton energy grows monotonically with the veloc- 
ity in the case of large a (see, e.g., a = 0.3 in Fig. |l|). In 
this case the soliton properties are qualitatively the same 
as in the limit of nearest-neighbor interactions (NNI's) 
for which Eq. (g) reduces to the Boussinesq equation. 
It is well known that this equation has a sech-shaped 
soliton solution w{z) = —1.5 (v^ — c^)/cosh {az), where 
a = •\/3(u^ — c^) is the inverse soliton width. As indi- 
cated above, the energy H '-^ [y'^ — c^Y^'^ of these solitons 
is monotonic function of the velocity, which means that 



there is only one soliton state for each given value of en- 
ergy or velocity. 



H 




FIG. 1. Energy of the pulse solitons versus velocity found 
numerically for different values of a and J (the value of J was 
chosen to get constant c — 1.515; see stars in Fig. bl): a = 0.3 
and J = 0.05 (dot-dashed line); a = 0.17 and J = 0.0172 
(long-dashed line); a = 0.1 and J = 0.0062 (dashed line); 
a = 0.05 and J = 0.0016 (fuU line). 
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FIG. 2. Shapes of pulse solitons found numerically for dif- 
ferent velocities at a = 0.17 and J = 0.0172. 

In the case of small a the soliton properties become 
much more interesting y,|7| . As it was recently shown [[7| , 
two branches of stable supersonic pulse solitons should be 
distinguished in this case: low-velocity and high- velocity 
solitons, separated by a gap with unstable soliton states. 

The solitons of the low- velocity branch are broad (they 
have a width much larger than l/s-|_), and can be de- 
scribed by Eq. (|5|) in the approximation (9^ — s^)w « 
— s^w. In this approximation the soliton solutions ex- 
ist in a finite interval of velocities, c < v < Vcr — 
\J['^(? — l)/3, and change their shape from the sech-form 
at II > c (see the case u = 1.6 in Fig. g) to the crest-form 
w(z) ~ exp(— a|z|/2) for i; -^ v^r (see the case v = 1.652 
in Fig. 0). Such crest solitons (or peakons) were first 



introduced in the theory of shallow water motion ||11|,|12[ . 
The solitons of the high-velocity branch are made up 
of two components: w = ws{z) + wl{z), where the short- 
range component 



wsiz) 
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is dominant in the center of the strain, while the long- 
range component wl(-z) — — 7 (s+''^^/12) exp(— s_|z|) is 
dominant in the tails (see the case u = 1.75 in Fig. ||). It 
should be stressed that this division of the soliton body 
into two components is not just a mathematical trick. 
Our present numerical simulations testify that the soli- 
tons of the high-velocity branch can be considered as 
bound states of the short-range and long-range compo- 
nents: they can be excited such that the relative distance 
between the components oscillates. However, such inter- 
nal soliton oscillations are highly damped and should not 
play an important part in the nonlinear dynamics of the 
system. The parameter 7 is determined by the equation 



7VA72 - 2{A + 2)7/3 + 1 = 3(A - 1) — v/1 - 27 (8) 

with A = K^/s'^. This equation, derived in Ref. |7|, has 
been there analyzed for small values of J and a, for which 
it has a unique solution at all values of velocity v. It has 
been shown that the interplay of the components ws{z) 
and u'l(z) results in this case into nonmonotonic depen- 
dence of soliton energy H on the velocity (see, e.g., the 
case a = 0.17 in Fig. ||), so that there is an energy in- 
terval where three soliton states with different velocities 
exist for each given value of energy. As is shown in Ref. 
0, the low- and high- velocity states (with dH/dv > 0) 
are stable while the intermediate state (with dH/dv < 0) 
is unstable. 
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FIG. 3. Three regions of the system parameters with qual- 
itatively difFerent properties of the pulse solitons. Full lines 
represent Eqs. IM) and (lid) . The points, marked off as circles 
and diamonds, were calculated numerically. Stars mark the 
parameters used in Fig. |l|. 



To sum up the foregoing, there is a demarcation line 
Ji(a) which separates the plane {a, J} into two regions 
(see Fig. y), namely: the Af- region (with a monotonic de- 
pendence of soliton energy on the velocity) at J < Ji{a) 
and the TV-region (with a nonmonotonic dependence of 
soliton energy on the velocity) at J > Ji{a). Our nu- 
merical calculations (see Fig. ||) validate the following 
estimation for Ji (a) : 



Ji(a) -0.23- 



(9) 



with ai ~ 0.25. The spectrum of stable soliton states 
is continuous and covers all supersonic velocities in the 
Af -region, while it has a gap (an interval of velocities 
with unstable soliton states) in the A'^-region. Emerging 
at J = Ji{a) this gap increases initially with growth of 
J. However, closer analytical examination of Eq. (0) 
shows that subsequently this gap starts to decrease and 
disappears again at J = J2{a) > Ji(a), where 



Ma) - o 
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with a2 — 0.16. Besides Eq. (hOf) we have found analyti- 
cal expressions for the soliton energy and impulse, all in 
a very good agreement with the numerical calculations. 
But due to lack of place we do not present these cum- 
bersome formulas in the present short paper. Instead, 
we just discuss below the results obtained for J > J2{a) 
with the intent to demonstrate that the soliton features 
in this region (lets call it Z-region) are qualitatively dif- 
ferent from those in the A'^-region. 
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FIG. 4. Shapes of pulse solitons which coexist at the same 
velocity v = 1.64, where a — 0.05 and J = 0.0016. These 
solitons are also indicated on the inset of Fig. hi 

Indeed, when J exceeds J2{a) there appears an interval 
of velocities in which Eq. (Q) has two real solutions. They 
correspond to two different types of two-component pulse 
solitons which coexist at the same velocity (6 and c in Fig. 



Q) . Accordingly, the dependence of the sohton energy on 
the velocity for J > J2(Qf) is not merely nonmonotonic 
but takes on a Z-shaped multivalued form (see, e.g., the 
cases a — 0.1 and a = 0.05 in Fig. |^). The possibility of 
such a dependence has been predicted in Rcf. using a 
variational approach. At that time, however, this predic- 
tion was met with disbelief and considered as an artifact 
of variational approach. But as we prove numerically in 
the present paper, the Z-region really exists. In this re- 
gion there is an interval of velocities where three soliton 
states of quite different shapes (see Fig. ||) and energies 
coexist at the same velocity. The soliton state with inter- 
mediate energy (b in Fig. 0) is always unstable. But the 
high-energy and low-energy solitons (a and c in Fig. |j) 
are usually (when dH/dv > 0) stable. The high-energy 
soliton on the low- velocity branch is broad and has a sin- 
gle component. But the low-energy soliton on the high- 
velocity branch and the soliton state with intermediate 
energy both consist of two components, short-range and 
long-range ones. The coexistence of two different types of 
stable pulse solitons at the same velocity causes new in- 
teresting phenomena, e.g., the synchronous propagation 
of two solitons with quite different widths (see Fig. 0). 



exists an interval of velocities where two types of stable 
pulse solitons coexist at each value of the velocity. The 
high-energy soliton is broad and has only a single com- 
ponent whereas the low-energy soliton consists of two 
components, short-range and long-range ones, and can 
be considered as a bound state of these components. It 
should be stressed that this phenomenon occurs even for 
LRI's of very sm,all intensity (hundreds times less than 
the intensity of NNI's, as is seen on Fig. 0) if only the 
radius of LRI's is large enough. It is important that the 
coexistence of two types of stable solitons does occur not 
only for the Kac-Baker LRI's discussed in this paper; 
on the contrary, it is rather a common phenomenon. In 
particular, we also have shown that it exists in anhar- 
monic chains with weakly screened Coulomb interactions 
between charged particles. 
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FIG. 5. Demonstration of the dynamical stability of the 
two stable pulse solitons (a and c in Fig. 0) which propa- 
gate with the same velocity v = 1.64. Here a — 0.05 and 
J = 0.0016. 

In conclusion, we show that the properties of pulse 
solitons in anharmonic chains with the long-range in- 
teratomic interactions are conveniently mapped onto the 
aJ-plane (see Fig. y). One can recognize in this plane 
three regions with qualitatively different properties of the 
pulse solitons. The M- and iV-regions were distinguished 
and discussed in Refs. HJ^ whereas the Z-region (that is 
the region of ultra-long-range interatomic interactions) is 
proven to exist in the present paper. In this region there 
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